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1 Introduction 



Symmetries and symmetry breaking serve as a basic principle in modern physics. 
For instance, it is generally appreciated now that all the four basic interactions in 
Nature can be formulated as theories with gauge invariance. However, there are 
still more to be understood, such as how gravity may be quantized, as quantum 
gravity is believed, for instance, to have dominated for a period (thought short) after 
big bang and determined the evolution of Universe. There seems a need for larger 
symmetries to set up a framework for quantum gravity. There has been argument 
that a fundamental theory of quantum gravity will not involve a space-time metric 
A theory depending not on a choice of space-time metric is invariant under 
diffeomorphism transformation, and is known as topological theory. This topological 
symmetry rules out any local dynamical degree of freedom, and the number of states 
of such a system could be as few as one - the vacuum - only. 

Topological quantum field theories are studied in the context of low dimensional 
topology 0. In particular, the topological quantum Yang-Mill theory is used to in- 
terpret the Donaldson theory which is regarded as a key to understanding the four 
dimensional topology, and thus as a key to understanding the geometry structure of 
the space-time [Q. In the three dimensions, the topological quantum Chern-Simons 
gauge theory is shown to provide a natural framework of studying the Jones poly- 
nomial of knot theory |Q, and to have important implications for two dimensional 
rational conformal field theories Q . 

Progress is also made in understanding quantum gravity. It is shown that (2 + 1) 
dimensional general relativity is equivalent to a topological field theory - the Chern- 
Simons gauge theory for the Poincare group, and thus, as a quantum field theory, is 
exactly soluble [^. Various topological models relating to gravity in different space- 
time dimensions have been suggested in the literature 00] [3. On the other hand, 
since any theory one uses to describe physics at low energy and/or low tempera- 
ture involves a spacetime metric and local dynamical degrees of freedom, there must 
be some means that connects a topological theory and its low energy/temperature 
limit. In particular, if a topological quantum field theory describes a phase of unbro- 
ken diffeomorphism of quantum gravity, a metric and local dynamics might arise as 
some form of symmetry breaking. Such a possible topological phase transition might 
happen in the early stage of Universe when it was extremely hot and small in size. 
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In this Paper, following a recent work |jTO|, we examine a system that exhibits a 
topological phase transition. The system we work with is a (2 + 1) dimensional one 
with a charged massive vector matter coupled to a U{1) gauge field. Both the vector 
and gauge fields are governed by a Chern-Simons term. While the Chern-Simons 
gauge field does not carry local dynamical degree of freedom (and is known as a 
topological field), the vector matter field carries them due to its mass term. Moreover, 
as the mass term involves a space-time metric, the massive vector field theory and thus 
the interacting theory are obviously topologically variant. Now, however, if somehow 
the vector field becomes massless, the local dynamical degrees of freedom, the space- 
time metric, and the difference between the gauge and matter fields disappear. And 
then the system , now with non-Abelian gauge symmetry, turns out to be topologically 
invariant. Namely, a topological phase transition would happen. The mass of vector 
field naturally serves as the order parameter characterizing the possible transition. 
To realize it, we heat up the system and see there indeed exists a critical temperature 
at which the effective vector mass vanishes. 

Breaking and restoration of various symmetries at high energy and/or finite tem- 
perature have been known for a long time in quantum field theories. Notable theoret- 
ical observations include the finite temperature restorations of SU{2) x f/(l) gauge 



symmetry for the weak electric interactions [ll||12l|13| and the deconfinement tran- 
sition in the quark-gluon plasma (for a review, see for instance fl^). Here we provide 
an example of phase transition concerning topological and gauge symmetries. The 
transition, recognized as of the zeroth order, is thus far as we know the most dis- 
continuous - even the free energy density discontinues at the critical temperature. 
Our conclusions are drawn based on perturbation expansion over the Chern-Simons 
coupling to the second order. We also observe that, with the thermal fiuctuation, 
the effective Chern-Simons coupling in the symmetry broken phase is just a slowly 
increasing function of temperature T with g{T = 0) = Qr and g{T = Tc) = 2gr at the 
second order, as we shall see below. Therefore the perturbation expansion and our 
conclusions should be reliable in the whole range of temperature T, as long as the 
renormalized zero temperature Chern-Simons coupling g,. is small. 

In the next section, we review the model in some details with emphasis on sym- 
metries. Section 3 is for deriving the partition function in path integral formalism. In 
section 4, we calculate one loop corrections to two- and three-point correlation func- 
tions. From these we obtain the critical temperature for the transition, and verify 
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the reliability of perturbation expansion. We then calculate in section 5 free energy 
density to two loops. From it, the order of topological phase transition is determined. 
Section 6 is devoted to summarizing our results. 



2 The Model 

The model of interests is a massive complex vector field 5^ minimally coupled to a 
gauge field who's dynamics is governed by a Chern-Simons term, 

/ = ^ ^ (eM.A5;(9, - iga,)B,^ + MB^B^ + t^.^a^d.a^) , (2.1) 

where the three space-time manifold Vl has Lorentzian signature. This theory was 
studied in a different context (anyon) |l5lll6|. This is an Abelian gauge field theory, 
violating parity and time reversal symmetries. Though the B^ field involves a Chern- 
Simons term as well, its nature is rather different from that of field, because of its 
non-vanishing mass term. Under a U{1) transformation parametrized by a(x, t), the 
Chern-Simons field varies as a gauge field 

a^(x, t) a^(x, t) + 9^a(x, t) , (2.2) 

and the massive B^ field transforms like a charged matter field 

i?^(x,t)-.e^"(^'*)i?^(x,t). (2.3) 

As a consequence of the U{1) symmetry, the Chern-Simons gauge field does not 
carry dynamical degree of freedom, but the massive vector field B^ does. To see this, 
we consider the equations of motion. For the Chern-Simons field a^, it is 

(-ixvxdyax = 3^ , (2.4) 

where the current = —ie^y\BlB\. From (^^ ), one obtains the current conservation 

^mJm = • (2.5) 



Moreover, from (|2.4|) , one can solve by expressing it in integrals over the current 
j^, with a use of gauge fixing to a^. This implies the Chern-Simons field is completely 
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determined by the matter field. On the other hand, the equation of motion for the 



massive vector field is 



e^^uxid, - iga,)Bx + iMB^ = . (2.6) 
Acting {d^ — iga^) on ( p.6|) , for M 7^ 0, one obtains 

(d^ - iga^)B^ = . (2.7) 



As is not an independent field, ( p.7|) ensures that one of the three components of 
B^ is eliminated in a covariant way. Another way to see the difference between and 
B^ fields is via the canonical approach, in which an independent degree of freedom 
must be accompanied by a non- vanishing canonical momentum conjugate, and vise 
versa. Let's write down the canonical momentum conjugate of the field variable 

TT^ = 7^ = TT^OMAflA • (2.8) 

Now one sees Qq is not an independent degree of freedom as ttq = 0. Moreover, the 
gauge transformation (|2.2|) implies one gauge degree of freedom of should be fixed. 
A simple choice is the axil gauge a2 = 0. Then, as ai's conjugate tti = = 
[see ( |2.^ )], ai is not an independent degree of freedom either. Therefore, the Chern- 
Simons field has no local dynamical degree of freedom at all. On the other hand, 
the canonical momentum conjugate of the massive -B^ field takes the same form as 
given in (|2.8| ). Therefore Bq is not an independent field variable. However, there is 
no gauge freedom associated to the massive B^ field, and Bi is the independent field 
variable and B2 the canonical momentum conjugate, or vise versa. 

On the other hand, when M = 0, the difference between and B^ disappears, 
except the latter is complex. Then, for convenience, we set gB^^ = Aj^ + iAj^ and 
gttfj, = y4°. Substituting these into (|2.1| ), we arrive at an action 



iiva Jn 

^2 



res = X ('^^^'^^" + le-'^AI.AlAlj , (2.9) 



where g = Arra. Carrying no local dynamical degree of freedom, this is known as a 
topological field theory. Now, the gauge invariance of the theory is not simply f/(l), 
but a non-Abelian one. The gauge group can be SU{2) or I SO {2, 1), with structure 
constant e°''"^. The variation of ^4° under a gauge transformation is 

= D^r» , (2.10) 
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where D^r" = d^r" + e"''"'A'l^T''. 

Moreover, independent of a choice of spacetime metric, (|2.9| ) is invariant under 
diffeomorphism transformations. Such a transformation can be generated by a vector 
on the three manifold Q, V^, via Lie derivative Cy- Under it, the Chern-Simons field 
An transforms as 

CvA; = D^iA^vn + VF^^ , (2.11) 

so that 

CvIcs[A] = 0, (2.12) 

provided vanishes on the boundaries of Q or the space-time manifold has no 
boundary 0||l3- However, subject to the flat connection condition F^^ = [_D^, D^] = 
0, diffeomorphism transformations completely fall into gauge transformations. In fact, 
the first term of ( |2.11D can be identified as the gauge transformation ( |2.1CI| ) with gauge 
parameter = A'^V'^, and the second term is proportional to the fiat connection 
condition. Therefore, on the constraint surface, a generator generates simultaneously 
both diffeomorphism and gauge transformation. Consequently, a gauge choice fixes 
both the gauge and diffeomorphism. As a topological quantum field theory, (|2.9|) is 
exactly soluble @]. 

Now, we see that the mass parameter M plays the role of order parameter. When 



M > 0, the system described by (|2.1| ) has U{1) symmetry only; however when M = 0, 
the system is invariant under non-Abelian gauge and topological transformations. In 
the following sections, we shall consider the quantization of the massive theory, i.e. 
the system in a symmetry broken phase, in a heat bath. We shall see then that the 
mass of the vector field becomes temperature dependent, M = M{T), and there 
exists a critical temperature Tc at which M(T) = 0, and a topological phase transition 
of the zeroth order happens. 

Another parameter of the theory is the Chern-Simons coupling g. It is dimen- 
sionless as the Chern-Simons interaction operator is marginal. This implies that the 
Chern-Simons coupling may need potentially non-trivial renormalization. However, 
due to the topological nature of the Chern-Simons term, g is not sensitive to changes 
of energy scale. In other words, the beta function of g vanishes identically, and g 
is not a running coupling constant. This fact makes a 'small' Chern-Simons cou- 
pling a perfect controlling parameter in perturbation expansion. On the other hand, 
however, when the thermal effect is taken into account, the concern is now that the 
Chern-Simons coupling may be temperature dependent, and if it went up rapidly with 
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temperature, perturbation would break down. We shall see below that fortunately it 
is not the case. The effective (temperature dependent) Chern-Simons coupling is just 
a slowly increasing function of temperature, and at the critical temperature it is only 
twice as large as the renormalized zero temperature one. 



3 Partition Function 

Now we attach the system to a heat bath with temperature T. As is known, finite 
temperature behavior of any theory is specified by the partition function 

Z = Tre-^^ ; (3.1) 

and the thermal expectations of physical observables 

< O > = lTr[Oe-^^] , (3.2) 

Zj 



where /5 = 1/T is the inverse temperature (fc^ = 1) 

As our system is relativistic, we like to work out its covariant formalism. It is 
clear from the discussion in the previous section, the Chern-Simons theory involves 
both the primary (first class) and secondary (second class) constraints. Therefore, to 
derive the partition function in the path integral formalism, special cares are needed. 



We shall outline the derivation in the zero temperature field theory [|Ty|, and then 
switch to the finite temperature case. We shall focus on the Chern-Simons field a^, 
i.e. consider the theory 

C^cs = -'^efiuxa^duax + a^j^ . (3.3) 

Generalization to the massive field will be then straightforward. With the momen- 
tum conjugates vr^ given in ( |2.8| ), it is easy to check that the theory ( |3.3| ) is subject 
to the following constraints: 



Xo = TTo , (3.4) 
1 

— ( 

2 

X = eijOittj+jo. (3.6) 



Xi = - -eijaj , (3.5) 



The last one is from 5Ccs/5ao = 0. Then the total Hamiltonian density is 

n = nc + Xoxo + kxi + Ax , (3.7) 
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where A's are the lagrange multiphers, and the canonical Hamiltonian density Tic is 



'He = a^TT^ - Ccs (3.8) 
= doVTo + ai{TCi - ^eijaj) + ao{eijdiaj + jo) - ajj • (3.9) 

From (|3.7|) and ( |3.9|) , one sees that do, dj and in H can be absorbed by the lagrange 
multipher Ao, Xi and A, respectively. Doing so, ao disappears, and so one can discard 
TTo by setting Ao = 0. The Hamiltonian density is now 

H. = XiXi + Ax - ttiji . (3.10) 

Now two constraints remain. The first, Xj, is of second class. And the second, 
seems to be of second class too. However, upon a linear combination with Xi, a 
redefined constraint, SjVTj + eijdiaj + jo, turns out to be of first class. As a matter 
of fact, this first class constraint is due to the gauge freedom of the Chern-Simons 
field. A nice approach to deal with a quantum system with first and second class 
constraints is to replace the Poison brackets with the Dirac brackets so that the 



second class constraints can be absorbed into the measure of the path integrals [21 



and the first class ones are taken care by using gauge fixing terms. For this purpose, 
let us choose a fixed time t. The constraint matrix element Cij at the time t is 

a,(x,y) = {Xi(x),X,(y)} = -2e,,5(x-y) ; (3.11) 

and the inverse 

Qi(x,y) = le,,5(x-y). (3.12) 
Then the basic (equal time) Dirac brackets (also called star brackets) are 

{ai(x),aj(y)}z) = ^eijS{^-y), (3.13) 
{7r,(x),7r,(y)}z5 = 2e,,5(x - y) , (3.14) 
{ai(x),7rj(y)}i5 = r]ij6{x-y). (3.15) 

And the first class constaint satisfies 

{x(x),x(y)}D = o. (3.16) 

Now, to fix the first class constraint x, one introduces a pair of ghosts, c and c, and 
their momentum conjugates, / and /, and the momentum conjugate of the lagrange 
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multiplier A, a. With these new fields, the system has the global nilpotent (off-shell) 
BRST symmetry |22|[|g. The BRST charge 

Q = ( d^x{cx-tfa) (3.17) 



generates the BRST transformations for all the fields via the Dirac brackets. Now, 
the partition function is 

Z,, = N J [Vij^y^ff , (3.18) 

with 

Ieff = J d?x (iiiT^i + Xa + cf + tf + aiji + {$, Q}d) , (3.19) 
where $ is a gauge fixing function, and the measure 

[P/i] = {det\x^,X,\y^^l[^{x^)^^a,V7r,VXVaVcVcVfVf . (3.20) 

i 

The determinant factor in the measure can be absorbed into the normalization factor 
A^. And the integrations over vTj eliminate the second class constraint Xi- Moreover, 
as a covariant gauge fixing for the later use, the gauge function is so chosen 

$ = ,c Qpa - d,a?j + fX (3.21) 

that 

{$, Q}d = -^cr^ - icdidiC + adiai - iff - A(jo + 2eijdia,) , (3.22) 

where p with dimension in the unit of mass parametrizes the covariant gauge fixing for 
f/(l) symmetry. Now, integrating out all the momenta (cr, /, / and VTj) and denoting 
the lagrange multiplier A as in the partition function, we obtain finally 



Zcs = N J Va^VcDcexp (i J d^xi^ef^^xafid^ax + a^^jf, + icd'^c + ^(c?^a^ 



(3.23) 

Above, the ghost fields c and c do not interact with other fields but only serve to 
cancel the non-physical, in fact all, degrees of freedom of The Chern-Simons field a^. 

The derivation of the partition function for the vector B^ field can be done sim- 
ilarly. The only difference is that there is no gauge freedom and so no first class 
constraint for the -B^ field, since it is massive. 

With these preparation, it is readily to work out the functional integral repre- 
sentation of partition function at finite temperature T. The trick is rather simple: 
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to replace the time variable t with the imaginary time ir via a Wick rotation, and 
to explain the final imaginary time as the inverse temperature (3 = 1/T. Then the 
partition function of system described by (|2.1| ) is 



Z = N I imUVa^VBlVBxVcVcexp I- f dr ( d^xC\ , (3.24) 

u X \ Jo J J 

with the Euclidean Lagrangian 

i Mi 1 

= --^tMuxBJ.id^ - iga^)Bx + —B*^B^ - -e^^xa^d^ax + (<9^c)((9^c) + —{df,a^f . 

(3.25) 

Being vector bosons or ghosts with ghost number ±1, all fields in ( p. 241 ) are subject 
to the periodic boundary condition such as 

a^{(3, x) = a^(0, x) and B^{f3, x) = B^{0, x) . (3.26) 



From the partition function (|3.24|) , it is easy to work out the finite temperature 



Feynman rules. The Chern-Simons propagator in the Landau gauge and the vertex 
are 

DIAp) = '-^. and rJ,, = ^6,.A. (3.27) 
And the vector propagator is 

W+M-' ■ ^ ^ ^ 



According to the periodic boundary condition ( |3.26|) , the third component of 



momentum, the frequency, takes discrete values, ps = 27mT for integer n's. Besides, 
each loop in a Feynman diagram carries an integration-summation TJ^nl 'i^p/(27r)^ 
over the internal momentum- frequency (p, 27rTn); and at each vertex, the momentum- 
frequency conservation is satisfied. 



4 Critical Temperature 

We now consider perturbation expansion of the theory. In this section, the two and 
three point correlation functions are calculated to the second order. For our purpose, 
we shall set the external momenta- frequencies zero. To start, we cite a useful equation 



10 



that maps the discrete summation TJ2^=-oo f{P^ = 2'KTn) into continuous integrals. 
By a contour integral on a complex plane, one obtains 

°° 1 fco rj^ roo+ie rl ^ 1 

Tjljin - 2.Tn) ^ - _|/(.)+/(-.„+/_^^__ _|;(„ + ;(_,„__ . 

(4.1) 

The expression is correct as long as the function /(ps) has no singularities along 
the real axis. One advantage of this equation is that it naturally separates the 
the temperature independent piece from the temperature dependent one. As an 
application, let's consider a typical integration-summation 

d'^p 1 



j(M.r) = ri:/^ (4.2) 



(27r)2p2 + ]Vf2 

which will appear in loop calculations. By using ( [4.1|) , the frequency sum is converted 
into the integrals 

J(MT)= I ^ r-^-d^r^ 1 1 

^ ' ' 7 (27r)3p2 + M2 J-oo+i. 27r 7 (27r)2 ;z2 + p2 + ]vf2 - 1 ■ 

The first term above is ultraviolet divergent and a regularization is necessary. If a 
naive cutoff A is introduced, the result is 

d^p 1 AM 



(4.4) 

/a (277)3 p2 + M2 27r2 An ^ ' 

On the other hand, if a gauge invariance reserved regularization such as the dimen- 
sional regularization is used, one ends up with the second term in (|4.4| ) only, though 
the integral is power-counting linear divergent in three dimensions. Namely, the two 
regularization schemes differ one another only by a A-dependent term. This sort of 
linearly cutoff dependent terms can be absorbed by re-definition of the zero temper- 
ature mass and coupling constant, as we shall see below, and therefore no physics 
should be affected by the regularization procedure(s) used. The second term in ( ^4.3| ) 
involves no divergent, thanks to the Bose-Einstein distribution function. This is a 
well-known feature of the finite temperature theory that all the divergences appear 
and thus can be taken care in the zero temperature field theories, and thermal effects 
are finite at any finite temperature T |23|[|^. The integrals on the complex z plane 
and on the real two-dimensional p space in the second term of (|4.3| ) are readily to 
perform, and we obtain 
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Now we consider one loop corrections to the correlation functions with vanishing 
external momenta. The inverse two-point function of 5^ field is defined as 

G,,-\p) = G%-\p)-^,.{p), (4.6) 

where Il^^{p) denoting the self-energy. Then the effective mass of the vector field 
can be defined as 

M{g,T)6,, = G^,-\p = 0) . (4.7) 

In general, with only 0(2) symmetry in the space, the effective mass M{g,T) of a 
vector field in the longitudinal direction is not necessarily the same with that in the 
transverse direction. For the massive vector field we are considering, it is the case. 

It is not difficult to check that the tadpole diagrams have not contributions to the 
self-energy. Then calculating 



S?J(P) =9'TT.J J^2'.^^Gl,{p + q)e,ruD'^^{q) , (4.8) 
we have 



S£)(P = 0) = ^A.J(M,T), (4.9) 



and the effective mass by using 

M{gr, T) = Mr + ^gl (m, + 2Tln(l - e-''^"")) + 0{g^) , (4.10) 

where denotes the renormalized zero temperature mass. For instance, = 
M — -^g'^A at one-loop in the regularization by a naive ultraviolet cutoff A, or 
Mr = M in the dimensional regularization. In the bracket in ( |4.10| ), the bare mass 



M has been replaced by the renormalized Mr, and the bare zero temperature Chern- 
Simons coupling g by the renormalized gr (see below), this replacement affects at 
most higher orders. The first term in the bracket of ( [4.10| ) is the radiative correction 
at zero temperature. Without a symmetry to restrict the renormalized vector mass 
to particular values, especially zero, the vector mass is a free parameter. Indeed, as 
we have seen here that, even one starts with a massless theory M = 0, the radiation 
correction at zero temperature generates a vector mass via the Coleman- Weinberg 
mechanism p5|. The second term in the bracket of ( [4.10| ) is obviously due to ex- 



changing energy with the heat reservoir. The low temperature limit T ^ is trivial, 
as M(g, T) (1 + -^g'^)Mr > 0. On the other hand, since M{g, T) is a monotonically 
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decreasing function of temperature, as T goes up, the thermal fluctuation tends to 
drive the effective mass to zero. Namely, there must exist a critical temperature Tc 
at which M{g, T) = 0, and a phase transition happens. Tc is readily to solve from 
( |4.10| ). We obtain at this order 

^-aMr/Tc ^ i_ ^-Mr/n ^ (4-11) 

with a = (l/2 + 37r/5f^). In a linearized form, as a good approximation when ^ M^, 

~ 3(^ + ^)Mr . (4.12) 

Now we see that the Chern-Simons interaction is responsible to the phase transition 
as it should be, and a stronger interaction causes a transition at a lower temperature, 
with the renormalized mass fixed. 

As is mentioned in section 2, the Chern-Simons coupling g receives only trivial 
correction at zero temperature so that the beta function of g identically vanishes. 
Alternatively said, it is insensitive to the change of energy scale. However, the thermal 
effect at finite temperature will affect the strength of the Chern-Simons coupling. 
Namely, in the finite temperature field theory, the coupling constant is a function 
of temperature, g = g{T). Then it is crucial to the perturbation expansion over the 
Chern-Simons coupling that g{T) must be reasonably small for T's in the range [0, T^]. 
The effective (finite temperature) coupling constant g{T) can be defined as 

r^.A(p = 0) = g{T)t^^.x , (4.13) 

where r^i,A(p) is the three-point function. Calculating the one loop diagram for the 
three vertex, we have 

Ca(p = 0) = ^(7%.aJ(M,T), (4.14) 
and then the effective coupling is 

9{T) = 9r{l~ ^9^.[l + 2^1n(l - 6"^^/^]) + 0(/) , (4.15) 

where the renormalized zero-temperature coupling gr = ^7(1 + 3^2^/ ) in the regular- 
ization by a naive cutoff, 01 g^ = g in the dimensional regularization. Again, in 
( |4.15|) , we have replaced the bare parameters g and M with the renormalized ones gr 
and Mri and this affects only higher orders. It shows that g{T) is a monotonically 
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slowly increasing function of the temperature T, with g(T = 0) = gr{l — -^gl) at zero 
temperature. On the other hand, at the critical temperature T = Tc, the effective 
coupling is 

g{T,) = 2gr . (4.16) 

This implies that in the region T < Tc, the perturbation expansion over the Chern- 
Simons coupling is reliable, only if the renormalized zero temperature coupling con- 
stant gr is small. 

Though the gauge fields at zero temperature have no dynamical mass due to 
the gauge symmetry, the quantum thermal fluctuations may endow them thermal 
masses. The electric and magnetic masses can be defined via the polarization tensor 
U^,ip,g,Mr,T) H: 

Meligr, Mr, T)5^o5,o + M,nag{gr, M„ T)5^j5,j = 'U^^ip = 0, gr, M„ T) . (4.17) 

Calculating the one- loop diagram for the polarization tensor at p = 0, we obtain 

ng)(p = 0,g,M,T) = ~g\. {j{M,T) + ^) , (4.18) 

and the thermal masses M.ei{gr, Mr, T) = A4mag{gr, Mr, T) = J^{gr, Mr, T) with 

M{gr,Mr,T) = -^glT\n{l - e"*^^/^) + 0{g') . (4.19) 

Above we have set the renormalized zero temperature masses of the gauge field to 
zero by using counter terms when it is necessary (for instance in the regularization 
by a large momentum cutoff), so that the gauge symmetry is respected. And we 
have replaced the bare zero temperature coupling g and vector mass M in (|4.19|) 
with the renormalized ones. The thermal masses of the gauge field are monotonically 
increasing function of T. At the critical temperature 

Min) = M,(l + ^g'^) + 0{g') . (4.20) 

OTT 

The electric and magnetic masses M.ei{T) and M.rnag{T) are known as the inverse 
screening length in the plasma [|I4[[|18[. Unlike those in the (3 + 1) dimensional QED 
and QCD where the magnetic screening is absent, our results here suggest that both 
the static electric and magnetic fields are screened by the plasma thermal excitations 
in the (2 + 1) dimensional massive system. And thus, at least to the one loop order, 
the plasma of thermal excitations acts like a superconductor. 
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5 Zeroth Order Transition 



We turn to calculation of free energy density for the system. As is known, the 
free energy is the single most important function in the thermodynamics. From 
it all other thermodynamic properties can be determined. In the present case, the 
free energy may tell us the nature of the topological phase transition. With the 
conventional Fourier transformation, we come to the momentum space. At the lowest 
order, performing the Gaussian integrals in (|3.24|) , we have the partion function 

Zo = [det(-e^,APA + M6^,)]-^[det{-e^,xPx + ^P^.Pu)]~^det{p^) . (5.1) 

The determinants are readily calculated. The one for the CS term, the second in the 
above, is 

[det(-e,.APA + -p,Pu)]-'^ = v(^nn(/^V)-' ■ (5.2) 

P n p 



The products are canceled by the determinant for ghost term, the last in ( |5.1| ), and 
the gauge parameter term ^/JJp contributes the zero-point energy. This verifies Chern- 
Simons gauge field carries no local dynamical degree of freedom. On the other hand, 
the determinant for the free massive field gives 

[det(-e,.,p, + M5,.)]-i = ^UWip" + M'r' • (5.3) 

Then the leading contribution to the free energy density, jFg = —IhiZq/^PV), is 

^o^^Ei:>"(^)+2r/|E,[,.^Ml-e-)] , (5.4) 

where u = \/p'^ + M'^. The second term is just the free energy of a gas of noninter- 
acting, massive bosons. Namely that the massive vector i?^ obeys the Bose-Einstein 
distribution. The factor "2" in the second term in ( |5.4| ) indicates two degrees of free- 
dom, carried by the complex field i?^, in the thermal equilibrium. Besides, like the 
gauge fixing parameter p, the vector mass M has an extra contribution to the zero 
point energy, as shown in the first term of ( ^.41) . 

Dropping the zero-point energy, the free energy density JFq can be re-written as 

^0 = --T'hA + 0{g') , (5.5) 

TT 1 
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where 

hn{x) = —— / dy ^— . (5.6) 

i \ n) Jo V2/ + X e V ^ "'"^ — 1 

The perturbation corrections to the free energy may be calculated by expanding 
the partition function (|3.24|) in the Chern-Simons coupling g and calculating the 



resulting Feynman diagrams. At the second order, it is 



InZo 



with the real (dashed) line standing for the (a^) propagator. Calculating the 
two-loop diagram, we obtain 

lnZ2 = -g^Mf3V {J{M, T)f , (5.7) 

and the correction to the free energy density 

^2 = J^,9lM^T^ + 21n(l - e-^^^/^)) ' + 0(/) . (5.8) 

Being positive-definite, J^2 increases monotonically with T. Physically, this implies 
the quantum thermal fluctuation tends to decrease the pressure -P(T) = — JF(T), 
contrary to the thermal behavior of the system in the free theory limit, as shown in 
( |5.5| ). It would be interesting to calculate the higher order corrections, for instance, 
the next order from the ring diagrams. To our main purpose of looking at the 
nature of the phase transition at large T ~ Tc, the result upto g^ is sufficient. 

At the critical Tc, as Tc ~ 3(| + p-)Mr. ^ Mr for weak CS interactions, we consider 
the leading term in both J-'q and JF2. It is readily to see 

J^oiTc) = -8l7r^h,{0)M^\ + O{h , (5.9) 

9r 9r 



^2{Tc) = \Ml^^+0{gl), (5.10) 



where /i4(0) = 2^(3) ~ 2.404. This shows that near the critical temperature Tc, 
JFq dominates. In other words, the pressure of the system to the second order is 
positive-definite as it approaches the critical temperature from the below, 

p{T ~ Tc) = -J^{T r^Tc)>0. (5.11) 
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Now let us look into the symmetry phase, where the mass of field vanishes. 
Due to the topological nature, the free energy must identically vanish. This can be 
readily verified within the covariant formalism used in this work. The symmetry 
phase of the system is described by the pure non-Abelian Chern-Simons theory ( |2.9| ). 
In a covariant gauge fixing, as usual one introduces ghost term 

d.c^D^c'' , (5.12) 

along with a gauge fixing term such as {d^A'^)^/{2p). Doing so, one fixes diffeo- 
morphism too. This ghost term includes the kinetic of ghost and gauge interaction 
between the ghost and Chern-Simons field. The Gaussian integral of the kinetic term 
of ghost with a = 1,2,3 cancel that of the Chern-Simons term for A^, and so 
JFq = 0. And then, order by order, the contributions to quantities like the free energy 
from the ghost loops cancel out completely those from the Chern-Simons field loops, 
at zero temperature as well as in the thermal ensemble. In particular. 



P(T) = . (5.13) 

This seems a natural consequence of lacking of dynamical degree of freedom in a 
topological theory (or in a topological phase of a system). 

Since the the pressure (or free energy) discontinues at the critical temperature Tc, 
we conclude that the topological phase transition is of the zeroth order. 



6 Summary 

In conclusion, we have looked into the thermal behavior of the theory with a charged 
massive vector matter coupled to Chern-Simons gauge field. To the second order 
in perturbation expansion, we have obtained a critical temperature Tc, at which the 
effective mass of vector field vanishes, and the system transfers from a symmetry 
broken phase to topological phase. The phase transition seems to be of the zeroth 
order, as the free energy of the system is discontinuous at Tc. We have seen that the 
effective Chern-Simons coupling g{T) is a slowly increasing function of temperature T 
in the symmetry broken phase, and therefore perturbation expansion over the Chern- 
Simons coupling should be reliable in the whole range of T, provided the renormalized 
zero temperature coupling is small. 
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The present investigation presents to our knowledge the first example that exhibits 
restoration of topological invariance at some critical temperature. At the phase tran- 
sition point, the local dynamical degrees of freedom and spacetime metric disappear. 
Meanwhile, the thermodynamical properties of the system subject to sudden changes. 
For instance, the pressure of the system promptly withdraws. And it is not difficult 
to check that (all) other thermodynamical quantities derivable from the free energy 
are discontinuous at the critical as well. This is a phenomenon that has not been 
observed in the Nature nor in theories, and deserves further studies. 

It is interesting to notice that the pure Chern-Simons theory with the non-Abelian 
gauge group IS0{2, 1) is equivalent to the (2 + 1) dimensional general relativity with 
a vanishing cosmological constant |Q. This is realized by identifying the vierbein 
and spin connection as the gauge field A^, and mapping the Einstein action to the 
Chern-Simons action (|2.9|) . Then, in this (2 + 1) dimensional theoretical laboratory, 
we have seen a possible realization of the topological invariance in quantum gravity. 
Namely that the thermal fiuctuation plays a key role to drive the dynamical system 
under consideration to a phase of pure gravity without a spacetime metric and local 
dynamical degree of freedom but with topological symmetry. 

The author thanks L. Dolan, C. Itoi, G. Semenoff and Y.S. Wu for discussions. 
This work was supported in part by the U.S. DOE under contract No. DE-FG05- 
85ER-40219. 
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